Generation of coherent terahertz pulses in Ruby at room temperature 
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We have shown that a coherently driven solid state medium can potentially produce strong con- 
trollable short pulses of THz radiation. The high efficiency of the technique is based on excitation of 
maximal TIfz coherence by applying resonant optical pulses to the medium. The excited coherence 
in the medium is connected to macroscopic polarization coupled to THz radiation. We have per- 
formed detailed simulations by solving the coupled density matrix and Maxwell equations. By using 
a simple V-type energy scheme for ruby, we have demonstrated that the energy of generated THz 
pulses ranges from hundreds of pico- Joules to nano- Joules at room temperature and micro- Joules 
at liquid helium temperature, with pulse durations from picoseconds to tens of nanoseconds. We 
have also suggested a coherent ruby source that lases on two optical wavelengths and simultaneously 
generates THz radiation. We discussed also possibilities of extension of the technique to different 
solid-state materials. 

PACS numbers: Gy.42.50 



I. INTRODUCTION 



The search for efficient, high-power, inexpensive, com- 



^— oact, and room-temperature methods of generation of co- 
lyRierent teraherz (THz) radiation is one of the main topics 
f^n modern optoelectronics and photonics ij. Its impor- 
^^ance is based on the fact that THz radiation has unique 
^^otential for a wide range of applications from diagnostics 
(-^ i different materials (including semiconductors, chemical 
Q^jpmpounds, biomolecules, and biotissues), imaging (for 
I 'medical and security purposes), to remote atmospheric 
Reusing and monitoring, astronomy, etc. 

In previous works we have focussed on a new approach 
the problem of generation of short coherent THz pulsed 
rTadiation by taking advantage of dramatic enhancement 
r^i the nonlinear response of a medium via maximal quan- 
'kPliim coherence |M El ; created in atomic and molecular 
*V&ases with different level configurations (double-A, V-A, 
cdiouble-V schemes) by coherent laser radiaiton. 

As is well-known, quantum coherence can dramati- 
cally change the optical properties of media. For exam- 
ple, lasing without population inversion (LWI) H and 
electromagnetically induced transparency (EIT) ^ have 
been demonstrated |3, ^| in coherently driven media. 
These unusual properties are employed for generation of 
electromagnetic radiation of different frequencies ranging 
from IR 1 1.1 T2j to UV ^1^] , and gamma-rays [ijj . En- 
hancement of coherent Raman scattering via maximal co- 
herence has been demonstrated experimentally |15| . and 



it has applications to enhanced real time spectroscopy 
(FAST CARS) 16]. Previously, the double A scheme 
with near-maximal Raman coherence was used for highly 
efficient conversion of blue to ultraviolet light in Pb vapor 
[l3l |. A new type of gas-phase optical parametric oscil- 
lator [Tll | was suggested for frequency downconversion to 
generate 1.88 /xm radiation in Pb vapor, and efficient in- 
frared upconversion was suggested to convert infrared 
light with a 100 /im wavelength to the visible. The first 
experimental demonstration supporting these ideas was 
performed in 4] where 5 /im IR radiaiton was generated 
in Rb atomic vapor. 

Many results in this area have been obtained in gases, 
but, recently, EIT has been generalized and extended to 
solids with a long-lived spin coherence a class of 

solid materials, namely, rare-earth and transition metal 
ion doped dielectrics have been suggested as very attrac- 
tive from the point of view of realization and applications 
of EIT. Several applications of EIT to improve perfor- 
mance of solid state lasers have been already proposed 
[TsL . LWI, first realized in a gas medium , recently 
has been demonstrated in solids as well [20l |. 

In this paper we extend the approach of generation of 
THz radiation via resonantly induced coherence to solid 
state media. The coherence at the corresponding transi- 
tion can be induced by coherent optical fields as shown in 
Fig.n Electronic levels of the THz transition are coupled 
by a pair of optical fields with Rabi frequencies Qi and 
to a common ground state forming a V system of energy 
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levels. Coherence Pcb induced by the optical fields at the 
THz transition gives rise to polarization provided that 
the corresponding dipole moment is not zero, which will 
radiate out a THz pulse with Rabi frequency f^a. This 
method differs from the one proposed in 5, 6] for gaseous 
media where the optically prepared THz coherence pro- 
duces no macroscopic polarization of gaseous media at 
a THz transition (this is the manifestation of symmetry 
with respect to inversion), the reason is that an electri- 
cally allowed dipole moment at the two-photon transition 
is zero since it is typically forbidden (although, it is im- 
portant to note here that the magnetic dipole moment 
can be non-zero even in gases, we discuss this opportunity 
later). For doped solids having sites with no inversion 
symmetry is common, such as the one occupied by Cr^+ 
(simplified energy level structure is shown in Fig.l^)) 
in ruby, the corresponding two-photon transition dipole 
moment is not zero, resulting in all three transitions in 
the V scheme being allowed. It is worth mentioning that 
using nonlinear optical mixing in a semiconductor het- 
erostructure for generation of few-cycle THz pulses was 
proposed in [23l |. 

In ruby closely spaced Ri and R2 lines arise from tran- 
sitions between the ground state of the Cr-^^ ion ('*A2) 
and its first excited state (^E). Each of these cubic-field 
states is split by the trigonal crystal field and spin-orbit 
coupling into a pair of Kramers doublets, 0.38 cm^^ apart 
in the ground state and 29.14 cm^^ apart in the ^E 
state. The width of the R lines is about 11 cm^^ at 
room temperature and 0.15 cm~^ at the temperature of 
liquid nitrogen. The transitions from the ^P^2 ground 
state to levels ^E are spin-forbidden, and have oscilla- 
tor strength ~ 7.5 • 10"'' (Ri, R2 line peak cross-section 
is ~ (1.2 — 1.4) • 10"^* cm^ at the temperature of liq- 
uid nitrogen and ^ 4 • 10"^" cm^ at room temperature). 
They are predominantly electric-dipole in nature, since 
the inversion symmetry of a cubic crystal field is bro- 
ken by crystal field trigonal distortions and odd-parity 
lattice vibrations. Namely, the odd-parity component of 
the crystal field at the Cr^+ ion mixes odd-parity states of 
high energy into the even parity d states between which 
R-transitions are observed and makes them weakly al- 
lowed. The magnetic-dipole contribution to the R-lines 
is estimated to be ~ 1 /30 of the strength of the elecric- 
dipole transitions [2^ . 

The 29 cm~^ transition in ruby has history which 
is as rich as the ruby laser itself. It was considered 
promising for realization of a quantum counter for far- 
infrared radiation (FIR) 27| on one hand, and as a sys- 
tem very convenient to study interaction processes of 



non-equilibrium phonons with two-level electronic sw- 
tems (so-called phonon spectrometer), on the other j2q |. 
There was also a proposal to make a far-infrared laser at 
this transition pumping via the R2 line, but due to the 
unfavorable ratio of relaxation rates (fast relaxation at 
the FIR and slow at the optical Ri transition), the laser 
could operate only at liquid helium temperature and the 
estimated gain did not exceed losses 29j. 

It is interesting that 29 cm~^ FIR was produced by 
non-linear mixing of two ruby laser beams, one emitting 
at Ri and another at R2, in LiNbOa 3Qj and in ZnTe 
crystals j^, with the efficiency 77 ~ 10~^. The method 
proposed in the present work would allow pulsed THz 
radiation to be easily produced from a dual-color ruby 
laser itself, operating at both Ri and R2 lines, thus sig- 
nificantly miniaturizing the system. 

The outline of the paper is the following: in the next 
section we describe a theoretical model used to calculate 
parameters of generated THz pulses. In section HI we 
consider possibilities of phase-matching for fields in dif- 
ferent regimes (CW and short pulses). In section IV we 
make estimates showing the perspective of the suggested 
scheme as well as discuss a proposed experiment with a 
standard femto-second laser setup. Next we show that 
the system we have considered is not limited, i.e. there 
are many systems to which our theory can be applied. 
Finally we compare the efficiency of our scheme with the 
methods of generation of THz radiation used nowadays 
and possible applications that are opened due to the high 
efficiency of the proposed technique. 



II. THEORETICAL MODEL 
A. Real system 

The ruby system shown in Fig.l^) includes eight 
states, four in the ground and four in the excited 
one. In our theoretical analysis we describe it using 
only three levels as in Fig.(^), since this simple model 
takes into account all essential physics and allows us to 
make estimates of expected THz radiation characteristics 
and required optical field parameters. 

For the proposed technique to work we need to make 
sure that an efficient V scheme can be realized. In the 
case of a single fs pulse, driving both optical transitions 
simultaneously, we will have only one polarization, for 
example, linear. For two pulses we can choose each pulse 
polarization separately. Selection rules for the right cir- 
cular and left circular polarizations are shown in Fig. 
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|24j . As one can see, neither for the right, nor for the 
left circular polarization a V scheme involving Ri and 
R2 transitions can be organized. For a linear polariza- 
tion perpendicular to the optic axis, which is an equal 
sum or difference of the right and left ones, there arc four 
V schemes possible, shown in Fig.ljJh)). 

Taking into account the relation between the coeffi- 
cients C+ and C~: C+C~* = -C~C+*, one can see that 
the schemes A and C cancel each other, since the prod- 
ucts of the matrix elements are opposite (±4C"'"C~*/9 for 
A and C, respectively). On the other hand, for schemes 
B and D the products are the same (-2%/2C+C^*/9). 
These two V schemes will interfere constructively and in- 
duce the THz coherence. 



B. Simplified model 

The interaction Hamiltonian for the system shown in 
Fig-©*) is given by 

Vi = -h [r2ie-"^''»> >< a\ + r22e^"^'="*|c >< a\ + h.c] 



h[n3e-''^-'''\c><b\+h.c.] 



(1) 



where f2,; = Pfiv^i/^ is the Rabi frequency of the respec- 
tive field, fiiy — ba, ca, cb and i — 1,2,3; pba and pca are 
the electric dipole matrix elements of the optical transi- 
tions b ^ a and c <-^- a, and Pcb of the THz transition 
c b, respectively; Wha, and Ucb are the frequencies 
of the electronic and THz transitions; £i is the amplitude 
of the respective electromagnetic field. 

The time-dependent density matrix equations are 



dp 
dt 



1 



H,p]--{Tp + pT), 



where V is the relaxation matrix, and H is the total 
Hamiltonian. 

To estimate the efficiency of THz radiation generation 
and required laser field intensities we first consider the 
system in the limit of ultrashort pulses, which are 
short compared to relaxation times of the system, and 
neglect propagation eflFects for the optical fields. In this 
simplified picture we assume that the system is driven 
by two resonant fields with equal time-dependent Rabi- 
frequencies Q,i = VI2 = ^{t) (assumed real). In the case 
of a V-scheme considered in this work the STIRAP tech- 
nique, proposed in 5, 6] for maximal molecular coherence 
excitation in a A energy system, is not applicable. In the 
case of the V scheme maximal coherence is excited by two 



fields of equal Rabi frequency and duration, acting simul- 
taneouly. The system is then described by the following 
set of equations 



dar 



dt 



{Pa - Pel - if^CTcfa, 



da. 



ha 



dt 



in [pa - Pb) - i^CFcb, 



dcFcb 

"df 



= {(^ba - ^ca) 



(3) 



(4) 



(5) 



dpa 

dt 



iVL {aba + <Jca) - in [<J*ba + <^*ca) > (6) 



dpb * s 



1 - Pa - Pb- 



(7) 



(8) 



Due to the symmetry of the model system pb — pc, (Jba — 
Uca and Ucb - real. Let us introduce new variables 

a = {(Jba + 0-ca)/2 = O-fca, (9) 



l3={a,b + cjl^)/2^Re{a,b) = (Tcb, (10) 

for which, taking into account lO and the following 
relation can be obtained 



- {pa + 213) = 0, 



(11) 



giving Pa + 2/3 = 1. For the new variables we arrive at 
the following system of equations 



drv 

— ^in- 4in(3, 

dt 



d£ 
dt 



2in (a* - a) , 



(12) 



(13) 



which can be further simplified if we substitute ^ — a — a 
d^ 



dt 



2in ~ 8in(3, 



dp 



-2in^. 



(14) 



(15) 
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The solution of this system is 



^ = ^ sin f2V2 C n{t')dt' ) , (16) 



/3 = = - ( 1 



cos 



2%/2 / n{t')dt' 



(17) 



It follows from Ea. p7|l that maximal coherence act — 
0.5 will be excited by a pair of pulses with area S = 
2\/2 D,{t)dt ~ TT each. Taking pulses of a Gaussian 
shape = nQexp{—t'^/2T^) with the area of the pulse 
S = Ay/nnQT, we will have the maximal coherence ex- 
cited when flQT = \/tt/4: sa 0.443. The excited coherence 
will give rise to the polarization which will radiate out a 
coherent THz pulse. It is worth noting that this mecha- 
nism is similar to free-induction decay (FID) |25| in pho- 
ton echo: excitation of a maximal coherence in a two-level 
system by a 7r/2 pulse, followed by emission of a coherent 
pulse. The difference is that typically in solids the FID 
decay is governed by the dephasing due to inhomogeneous 
broadening of the transition leading to decay of the corre- 
sponding coherence during the time ~ Wj^h- case 
of ruby the THz transition is homogeneously broadened 
even at liquid helium temperature, so the decay of the 
THz coherence is determined by the homogeneous width 
of the corresponding transition. 

For two long pulses, which can be considered as CW (in 
the sense that optical and THz coherences approach the 
steady-state), with approximately equal Rabi frequencies 
ill « ri2, resonant with the corresponding optical transi- 
tions, we get the THz coherence 



2^2^ 



* exp 



7+(|Oi|2 + |n2|2)/^^i, 



7cfc(7+(l^^iP + |f^2P)/7cfc) 



(18) 



where ^ba = lea — 7 and icb are the optical and THz 
coherence decay rates, respectively. In this case the max- 
imal THz coherence \acb\ = 1 will be excited provided 
that (|f2ip + |il2p) ^ 77cb, which is a usual steady-state 
EIT threshold condition. 

Let us now estimate the peak intensity and energy of 
the optical pulses necessary to excite required coherence. 
The peak intensity of the pulse is expressed in terms of 
the peak Rabi-frequency as 

_ 2nncnl 

"fnXaabs 

where 7 is the width of an optical transition, n is the 
refractive index at the optical wavelength, A is the opti- 
cal wavelength, and Uabs is the corresponding absorption 



cross-section. The energy of a single pulse can be calcu- 
lated as 



E. 



opt 



2'K^/T:ficA, 



opt ^2 



inXaabs 



(20) 



where Aopt is the laser beam cross-section. 



III. PROPAGATION: PHASE-MATCHING 
GEOMETRY FOR FIELDS 

A complete self-consistent calculation also includes 
field propagation equations 



dz 



dt 



(21) 



where index a — 1,2,3 indicates all fields and correspond- 
ing polarizations, rja — 2'irujaN^^/nach is the corre- 
sponding coupling constant, uja are the frequencies of the 
optical and THz fields, N is the density of the medium, 
c is the speed of light in vacuum, Hq, is the correspond- 
ing refractive index, and losses for the field during 
propagation in the crystal due to scattering, diffraction, 
or non-resonant absorption. For optical fields these losses 
are usually small, but for the THz field in free space the 
diffraction losses given by K3 — Xs/D^ should be taken 
into account. To avoid diffraction losses the crystal can 
be placed in a waveguide for THz radiation, then the dis- 
tribution of the field mode should be taken into account. 

Propagation effects are important for nonlinear inter- 
actions such as wave mixing. Let us consider two regimes: 
CW, when the spectral widths of pulses is smaller than 
the splitting between the levels of the THz transition, and 
pulsed one, when the spectral width of pulses exceeds the 
splitting between the levels. 

In the CW regime the field fia at the output of the crys- 
tal is given by (assuming that the THz pulse propagates 
in z direction) 



n 



3 = «e 



-K3L 



g tdk + K3 



(22) 



sin{SkL/2) 
{5kL/2) ' 



if losses are neglected. Here L is the length of the crystal, 
and the prefactor describes phase-matching (Sk = k2z — 
kiz — ks). The direction of THz propagation is given by 
the phase-matching condition (see Fig. 3) 



k2 - ki = fcs 



(23) 
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To satisfy the phase-matching condition 6k = for the 
THz pulse to propagate in z direction (along the crystal) 
the angle between the optical fields has to be adjusted 
according to the equantion 



k\ + — 2kik2 cos (j) = k^, cos 



2fcifc2 



(24) 



Rewriting a condition for optical and THz frequencies 
0^2 — wi = t03 in the form of ks = (fe2 — ki)n3/ni, we 
obtain 



(25) 



where Ai(3) is the wavelength of the optical (THz) field. 
For the case of ruby the refractive index for THz radia- 
tion is twice larger than the one for optics, « 2ni^2 
[29I I , then the angle between directions of propagation for 
optical fields = 3- 10"'^. The angle between the ki(^2) 
vectors and the THz field propagation direction is 7r/3 in 



this case. This provides an interesting opportunity to de- 
sign a dual-wavelength ruby laser working simultaneously 
at i?i and R2 lines that can also generate THz radiation, 
as is shown in Fig|31;. The crystal will simultaneously 
serve as a THz waveguide in this case. 

Now we turn to a short pulse excitation regime. For 
many femtosecond laser experimental setups it would be 
an interesting extension to include an option of THz gen- 
eration. The phase-matching in this case can be obtained 
by excitation of THz coherence with two optical beams 
with the same frequencies but different directions of prop- 
agation. It is instructive to show how it works in this 
case. Let us consider the problem in the simplest case of 
lowest order of perturbation in coherent pumping fields. 
The density matrix equations look like 



cr'ab = -iflie 



in 
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(26) 



(27) 



J 



n2{t)e 



ikxX--\-ikzZ\ 



A solution of this set of equations is given by 



ini[t)e 



"ik^ cc ik ^ 2^ 



n2it)e- 



(28) 



a^b = 2 [ dt' I dt" {ni{t')ni{t") + n2{t')n2{t" 



ni{t')n2it") cos2k,z + VLi{t")VL2{t')cos2k^z) 



(29) 



r 



One can see that the phase-matching between the optical 
fields and the THz field is achieved by appropriate tilting 
of the optical beams at the angle (p « 2kz/kx = k^/ki = 
Aina/Asni; in particularly, for ruby (j) « 10"'^. 



IV. ESTIMATION OF ENERGY OF 
GENERATED THZ FIELD 

First we analize the optimal conditions at which the 
maximal THz coherence can be induced in ruby by the 
optical fields at both room and liquid helium tempera- 



tures. For room temperature ruby 7 — ^cb = 5.5 cm ^ 
■3^. n = 1.76, A = 694.3 nm, and aabs = 4 • 10~2° cm^ 
|26| . In the case of ultrashort driving pulses (100 fs-1 
ps) the shortest THz pulse that can be generated in this 
system has a duration ^ 7^^ ^ 1 ps. With the optical 
pulses of the same duration t = 1 ps the Rabi frequency 
rio = 4.43 • 10^^ s~^ is required to produce the maximal 
coherence, corresponding to Ipeak ^ 10^^ W/cm^, which 
is still below a damage threshold of ruby for such duration 
of the pulses 33] , but we make a conservative estimate for 



lower intensities /„ 



100 GW/cm , certainly below 
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the threshold, corresponding to the peak Rabi frequency 
Oq = 2 • 10^^ s~^. For long driving pulses (10 ps -1 ns) 
to excite the maximal coherence \acb\ = 1 would require 
^1.2 ^ 77cfc resulting in the Rabi frequency Qi^2 ~ 10"^^ 
s~^, corresponding to intensities exceeding the ruby dam- 
age threshold for pulses of such duration. For long pulses 
(1-10 ns) the damage threshold intensity is in the range 
Ipeak -- 20 - 30 GW/cm^ '34], which corresponds to the 
Rabi frequency fli.2 ^ 10^^ s^^. This value is used in the 
estimates below. 

At low temperatures (1.8 — 4.2 K) the THz coherence 
lifetime is significantly larger, ~ 500 ps 1^3, which 
means that with ultrashort pumping pulses (in this case 
of 10 ps-1 ns duration), the maximal coherence will be 
excited with only fig ~ 10^° s~^, requiring Ipeak ~ 100 
MW/cm^ pulse intensity. For long pulses (1-10 ns) the 
maximal coherence will be excited when flf 2 11 cb, and 
since at low temperatures 7 '-^ 0.05 cm"-"^ |3^, it results 
in 0.1^2 ~ 5 • 10^ s"\ requiring Ipeak ~ 30 MW/cm^. 

The next question is the energy of the generated THz 
pulse, which can be estimated from the propagation equa- 
tion 



dz 



123 dfl3 _ 27rw3«iV^^j 



dt 



n^ch 



(30) 



where N is the density of Cr'^+ ions, huj^ = 29 
cm~^ (0.87 THz), the dipole moment of the far- 
infrared transition can be calculated using measured ex- 
perimentally low-temperature (LT) parameters /i^^ — 

Icb^'^^^s'^abs THz/'^'^^3' ^^'^ nou-resouant ab- 

sorption coefficient of the sapphire host at the THz wave- 
length. 

Assuming that the crystal sample is side-pumped ho- 



mogeneously, so that the THz coherence acb does not 
depend on the propagation coordinate z, and phase- 
matching conditions are satisfied we arrive at the equa- 
tion 



dz c dt 



irjUcbit) - K,^3, 



(31) 



where 77 = N'f^^cr^^^ thz^ so that at the end of the crys- 
tal the Rabi- frequency of the THz pulse is 



n^it, L)^i7jj aeb [t - ^^n^j e-<'-'')dz', 

(32) 

where Zq = max (0, L — tc/ns). For the ultrashort pump- 
ing pulses we can model acb{t) = cr^"^ exp {—'fcbt) at 
t > and acb = at t < for an estimate. Ea. (|32(l gives 



Icb-nzjc ■ 
for t < Lns/c and 



Q3 — irju 



max 
cb 



Icbns/c^ 



for t > Ln^/c. 

Eas. (|33|l - H34() then allow one to calculate the energy of 
the THz pulse 



THz 



Icb "^^^abs^THz Jo 



\^:^\'dt 



(35) 



maa:\2 
cb ) 



2K-icb{K + lcbn3lc) 2K-feb{lcbn3/c- k) 7^(, ((7cb"3/c)^ - 



where A is the crystal cross-section along THz pulse prop- 
agation direction. For two pulses of 1 ps duration with 
Ipeak ~ 100 GW/cm^ the maximal induced coherence is 
^rnax ^ o.21. Taking iV = 1.6 • IQis cm^^ for 0.05% 
doped ruby, and known from low-temperature measure- 
ments 7^^ = 2-109 s'S a^^^^, = 3 • lO'^^ cm^ H, 
and room temperature parameters jcb — 10^^ s~^ and 



K = 0.4 — 0.5 cm~^ [23, 713 = 3.5 j^^, and considering the 
crystal of 1 cmxl cmxO.l cm size with the A—1 cm x 
0.1 cm THz emitting cross-section and L = 1 cm length, 
we arrive at Ethz ~ 630 pJ. The energy can be higher for 
higher optical peak intensities, approaching several nJ for 
Ipeak ~ 10^^ W/cm^. For low temperatures (1.8-4.2 K) 
the parameters are: jcb — 2-10^ s~^, (Jabs,THz = 3-10^^9 
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cm^, K = 0.01 cm~^ and the maximal coherence 

fcfc ~ 0.5 can be easily excited. Due to stronger absorp- 
tion at low temperatures {aabs ^ 10~^^ cm^ ,26]), half 
of the optical intensity will be absorbed at « 0.5 mm for 
0.05% Cr'^^ density, therefore, for this estimate we take a 
crystal with dimensions 1 cmxO.l cmx0.05 cm size with 
the THz emitting cross-section A — 0.1 cmx 0.05 cm and 
L = 1 cm length which results in Ethz 7.5 /zJ. 

Let us also make an estimate for fs pulse duration. 
Since powerful fs Ti:Sapphire lasers are readily available 
in many laboratories nowadays this is of interest. In this 

I 



where from Eq.l^HIl ^S""^ = 

2n*2ni/ {-fjcb + Pil'' + W) and G 
6|f^i|V(7+(|f^i|' + |f^2n/7c6)- In the limit 

^^1 2 — 77c6 the maximal coherence is excited CTcfc = 1 , 
which would require, though, experimentally unfeasible 
Rabi frequency f^i 2 ^ 10^^ s^^ at room temperature, 
resulting in Ipeak ^ 10^^ W/cm^. For such long pulses 
(10 ps -1 ns) the damage threshold intensity for ruby is 
jpeak ^ 20 — 30 GW/cm^, which corresponds to the Rabi 
frequency rii_2 ~ 10^^ s^^. For this Rabi frequency and 
room temperature decay rates 7 ~ jcb = 10^^ s""'^, the 
ratio ilf 2/77cfc ~ 10~^ and the resulting THz coherence 
IfJcbl ^ 10^^ for a 10 ps pulse. The THz pulse energy will 
be Ethz ~ 300 pJ. At low temperatures the maximal 
coherence acb = 1 can be excited by long pulses (1-10 
ns), which would give energy of a THz pulse Ethz ^ 8 
/iJ. 

Finally we estimate the total energy of optical pulses 
using Ea. H20|l and overall conversion efhciency to THz 
radiation. For room temperature 1 ps pulses pumping 
the side of the crystal with Aopt = 1 cmx 0.1 cm size, 
the energy in one pulse according to Ea. lf^ is Eopt ~ 
27 mJ, for two pulses, 54 mJ, respectively. We have to 
take into account that only (1 — exp (— (TahsA^iopt)) of 
the incident energy is absorbed, where Lgpt is the size of 
the crystal along which the optical pulses propagate. For 
room temperature aabs = 4 • lO^'^" cm'^, N = 1.6 • 10^^ 
cm~'^ and Lopt = 1 cm, the absorbed fraction is 0.47, 
so about twice the energy calculated above is required, 



case the coherence can be excited by just one pulse of 100 
fs duration, since its spectral width is larger than the THz 
transition splitting. Although, to satisfy phase-matching 
conditions, it is necessary to split it in two beams hiting 
the crystal at slightly different angles. Keeping the same 
peak intensity Ipeak ~ 100 GW/cm^ we obtain a smaller 
pulse area S — 0.14, giving the coherence ct^S""^ ~ 2.5 • 
10"'^ and resulting in THz pulse energy of « 100 fJ at 
room temperature. 

For the CW pumping case we have, in analogy with 
the short pulse pumping. 



(36) 

I 

resulting in ^ 110 mJ total energy. For a pair of 100 
fs pulses the same reasoning leads to the total energy of 
~ 11 mJ (2.7 mJ in one pulse from Ea. (|20|l ). For long 
pulses the total required energy is 300 mJ. 

The expected values of the energy of THz pulses, the 
peak amplitude of the THz field if focused to a spot of 300 
/j,m size, required pumping radiation energy and duration 
and efficiency of optical to THz energy conversion for 
room temperature ruby are summurized in Tabled 



At low temperatures as was already discussed above 
55% of optical intensity is absorbed at Lopt = 0.05 cm. 
Taking Aopt — 1 cmx 0.1 cm we obtain from Ea. (|20|l for 
short pulses Eopt ~ 13.5 mJ in one pulse, giving total 
required energy of 54 mJ. For long pulses the estimate 
gives Eopt = 67.5 mJ in one pulse and 270 mJ total energy 
for long pulses. 

The parameters of THz pulses: energy, peak amplitude 
if focused to a spot of 300 /im size, and required optical 
pulses energy and duration at low temperatures are sum- 
murized in Table II. 



nW3^A^^7i;,^CT^ZrF/. 9/ 1 e-2KL -{K+Gn^/c)L 

TP 'CO abs,lHz /max\'i I ^ ^ ■ ^ 

t^THz - {CFcb ) 



ns ' ' \2nG{n + Gn:ilc) ^KGiGn^jc-K) g ((Gna/c)' 
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TABLE I: Estimates of the THz pulse energy, peak THz field amplitude, required optical pulse parameters and THz radiation 
generation efficiency in ruby at room temperature 



ultrashort 


pulses (100 fs-1 ps) 


long pulses (10 ps-100 ps) 


100 fs 


1 ps 


10 ps 


acb ~ 2.5 ■ 10'^ 


acb = 0.21 


acb ~ 10~^ 


Ethz ~ 100 fJ 


Ethz ~ 630 pJ 


£thz ~ 300 pJ 


Sthz ~ 300 V/cm 


Sthz ~ 23 kV/cm 


£thz ~ 5 kV/cm 


^0 = 2 ■ 10'' s"' 


^0 = 2 • 10'' s~' 


fio ~ 10" s"' 


Ipeak = 100 GW/cm^ 


Ipeak = 100 GW/cm^ 


Ipcak = 20 - 30 GW/cm^ 


Eopt = 110 mJ 


Eopt = 11 mJ 


Eopt ~ 300 mJ 


"77 ~ 6- 10"^ 


77 ~ 10~" 


ri ~ 10"* 



"Efficiency is estimates as a ratio of the THz pulse energy to 
the total energy of optical pulses absorbed in the crystal r] = 
^THz/'^Eopt- The ultimate efficiency of the method is given by 
the ratio of the THz radiation frequency to the frequency of the 
optical transitions 29 cm-'/14420 cm-'=2 • 10~^. 



TABLE II: Estimates of the THz pulse energy, peak THz field 
amplitude, required optical pulse parameters and efficiency of 
THz radiation generation in ruby at low temperatures 



ultrashort pulses (10 ps-1 ns) 


long pulses (1 ns-10 ns) 


acb = 0.5 


(Jcb = 1 


Ethz ~ 7.5 /iJ 


Ethz ~ 8 /xJ 


Ethz ~ 230 kV/cm 


Sthz ~ 170 kV/cm 


= 10'" S-' 


f7o ~ 5 ■ 10^ S-' 


Ipcak = 100 MW/cm^ 


Ipcak = 30 MW/cm^ 


Eopt = 54 mJ 


Eopt ~ 270 mJ 


7J ~ 10-* 


?7~3-10-5 



V. DISCUSSION 
A. Application to other solid materials 

It is important to note that the proposed technique can 
be appUed to other solid materials with suitable transi- 
tions in the THz range. Similar schemes can be found in 
other ions as well, since for a typical rare-earth or transi- 
tion metal ion the level structure is complex, having many 
crystal field and spin-orbit split components, separated 
by tens to hundreds of wave numbers. For example, this 
method is readily applicable to Cr'^+:BeAl204 (alexan- 
drite), having R-lines splittin g of 41 cm^^ (1.23 THz) 
in a non-inversion Cr"^+ site |3^, and R-lines absorp- 
tion/emission cross-section ten times larger than ruby 
(3 • 10"^^ cm^ at room temperature). It allows one to 
decrease peak optical intensities to induce maximal THz 
coherence and results in higher convertion efficiency. 



The technique can also be applied to rare-earth ion 
doped materials, for example Pr^+ doped hosts, such as 
LaFa YzSiOs [io], CaFa puj with the ground and 
excited state Stark splittings being in a broad range 17 
- 100 cm^^ and higher, the same splittings are found in 
Nd3+ doped crystals, such as YAG 'i^ and YVO4 0, 
and a number of others, for example Tm:YAG |4J| and 
Er:YAG 45] and Er:YLiF4 

Also it is important to note that the technique can be 
applied to gases as well, but the coupling to THz radi- 
ation in this case occurs by a magnetic-dipole moment 
(macroscopic polarization due to an electric-dipole mo- 
ment is zero because of symmetry reasons). Typically 
the magnetic dipole moment is smaller than the electric 
dipole one, so efficiency is smaller, but for higher den- 
sities and pressures this approach can still produce THz 
radiation. The dipole moments for optical transitions can 
be stronger than for ruby, so the coherence can be larger. 
This possibility is an extension of the method proposed 
in 0,13, although its analysis is beyond of the scope of 
the present paper. 

The next important remark is related to the case when 
the THz transition of interest belongs to the ground state. 
Then there are a lot of new opportunities to use more 
sophisticated methods of coherence preparation. In par- 
ticular, the method of pulsed coherence production via 
Stimulated Raman adiabatic passage (STIRAP) can 
be used to control the duration of the THz pulse similarly 
to 13]. Generally, all methods for coherence preparation 
via femtosecond pulse shaping developed to improve the 
sensitivity of coherent Raman scattering in partic- 
ular, frequency chirping and fractional STIRAP [T5| . 
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matched pulses 0|, etc. can be applied. 



B. Comparison with current methods of short THz 
pulse generation 

It is interesting to compare the method of THz genera- 
tion proposed in recent works and in this paper with 
currently available ones. There are a variety of methods 
already considered successful for generation of short THz 
pulses, so it is worth to compare the above estimates with 
the currently achieved parameters. 

At the moment the most impressive results in gener- 
ation of short THz pulses were achieved using electron 
beam based sources^ free-electron lasers (FEL) l48j 
and synchrotrons 49, "5^. For FELs, typically THz pulses 
of about 1 ps duration with 1-40 ns distance between mi- 
cropulses are generated, grouped into a few-^s trains. En- 
ergies of the micropulses are about 1-50 mJ (the param- 
eters of FELBE laser, Rossendorf, Germany). Recently, 
substantially higher powers of coherent broadband THz 
pulses, produced by synchrotron emission, were obtained 
from the electron beamline 49, 50] . Hundreds of fs-short 
half-cycle THz pulses were generated with energies up to 
100 AiJ. 

THz pulses with ns durations were achieved in THz 
semiconductor lasers (quantum cascade lasers, p-Ge and 
n-Si lasers [53|-[5^), limited by the need of cryogenic 
cooling and in the case of p-Ge lasers by high electric 
and magnetic fields required for laser operation. In p-Ge 
lasers 52J- 54] 20 ps pulse durations were obtained in a 
mode-locked regime with peak power up to several Watts 
in a few-/is train of pulses. 

The obvious advantages of the optical crystal based 
THz emitters discussed in this work over FEL and syn- 
chrotron sources are compactness (typical optical solid- 
state materials are significantly smaller), cost, and ease of 
handling. Compared to THz QC lasers the crystals are 
much easier to grow and handle, they also can provide 
better coupling between the THz mode and the generat- 
ing material due to a larger size of crystals. 

Other methods of short THz pulse generation are 
based on interaction of different materials with ultra- 
short laser pulses. The most popular ones, giving subpi- 
cosecond THz pulses, are the photocurrent method using 
the Auston-switch technique [Sg and the optical rectifi- 
cation technique [s^ . Since the proposed crystal-based 
THz sources are closer in their characteristics to these 
methods, we will make a more detailed comparison. Re- 
cently THz sources based on amlifier-laser systems such 



as Ti:Sapphire, utilizing different methods of THz ra- 
diation generation such as optical rectification in non- 
linear crystals and laser-produced plasma and photocur- 
rent THz generation in semiconductor antennas, became 
popular. The transition metal and rare-earth ion doped 
crystals proposed in this work as THz emitters are suit- 
able for amplifier-laser-based THz systems working in 
a single-shot regime 58], since the typically long (hun- 
dreds fis - milliseconds) population decay time of excited 
states of these ions limits the repetion rate. Crystals 
are also better in terms of saturation, able to withstand 
high laser pulse fluence and thus very compatible with 
low-repetiotion rate table-top laser systems producing 
pulses with energies in the 10-100 mJ range with ~ 10 Hz 
repetion rates. Fig^ compares the predicted in f5^ per- 
formance of different amplifier-laser-based THz sources 
and ruby at room temperature, at optical pulses energy 
of ~ 10 mJ. It shows that ruby is expected to perform as 
good as the biased GaAs antenna with 1 kV/cm dc bias 
field and the plasma THz source with an external bias 
field, slightly yielding them in conversion efficiency. 

FigEl compares the predicted in [E^ performance of 
different amplifier-laser-based THz sources at room tem- 
perature and ruby at liquid helium one, at optical pulses 
energy of ~ 100 mJ. The ruby source at liquid helium 
temperature can compete with the the plasma source uti- 
lizing the fundumental and second harmonic optical fields 
in terms of THz pulse energy and conversion efficiency, 
and yields it in peak THz field amplitude. 



C. Advantages of the technique 

New methods of THz generation proposed in 0, |^ |^ 
for atomic and molecular gases and in the present work 
for doped crystals potentially have very high efficiency, 
they offer the possibility of THz pulses with controlled 
durations including femtosecond region and the possi- 
bility of generation of pulses with high energy, already 
comparable with synchrotron or free-electron laser based 
THz sources. Implementation of the method, therefore, 
will open exciting opportunities in many fields of THz 
applications such as the time-resolved THz spectroscopy. 
Using intense short pulses of THz radiation generated by 
the proposed method will significantly increase the per- 
formance of spectroscopic measurements, allowing max- 
imal (among currently existing methods) temporal reso- 
lution. Extremely high efficiency of the proposed tech- 
nique opens a way for orders of magnitude increase in 
sensitivity of the spectroscopic methods, THz tomogra- 
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phy, non-destructive quality control, medical diagnostics 
and biomaterial characterization. 

Another advantage of the proposed scheme is a better 
control over the phase of the generated THz radiation, 
resulting in lower (pulse duration x bandwidth) product 
compared to the non-coherent techniques. Resulting THz 
pulses have narrower bandwidth, thus offering an oppor- 
tunity of better resolution measurements. Coherence of 
the generated THz pulses combined with their high en- 
ergy opens the way to various non-linear phenomena in 
the THz range, as well as to coherent phenomena simi- 
lar to the ones observed with lasers in the visible range, 
such as Rabi oscillations, coherent transients, etc. The 
proposed methods have, therefore, high potential for in- 
vestigation of nonlinear effects in chemical and biological 
objects, in medical applications, which is a substantial, 
practically unexplored field of research. 

VI. CONCLUSION 

In the present work we proposed a new technique 
for generation of short intense THz pulses in coherently 
driven doped optical crystals. The method is based on 
excitation of maximal THz coherence by a pair of reso- 
nant optical pulses, resulting in polarization build-up in 
the medium, which will radiate a THz pulse. As an ex- 
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FIG. 2: a) Selection rules and transition matrix elements for 
Ri and R2 optical transitions in ruby for right and left cir- 
cular polarizations; b) Four V schemes possible for linearly 
polarized fields. 




FIG. 3: a) CW fields 1 and 2 induce coherence between levels 
h and c in ruby; b) For short pulses fci = fc2; c) A three mode 
ruby laser that generates optical radiation at _Ri and R2 lines 
and THz radiation simultaneouly with the ruby crystal serving 
as a THz waveguide. 
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FIG. 4: Predicted in [53 performance of several amplifier-laser-based THz sources and ruby at room temperature with optical 
pulses energy ~ 10 mJ: a) peak THz field amplitude; b) THz pulse energy; c) conversion efficiency. 
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FIG. 5: Predicted in ^3 performance of several amplifier-laser-based THz sources at room temperature and ruby at liquid 
helium one with optical pulses energy ~ 100 mJ: a) peak THz field amplitude; b) THz pulse energy; c) conversion efficiency. 



